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ABSTRACT 

The detection of mixed oscillation modes offers a unique insight into the internal structure of 
core helium burning (CHeB) stars. The stellar structure during CHeB is very uncertain be¬ 
cause the growth of the convective core, and/or the development of a semiconvection zone, 
is critically dependent on the treatment of convective boundaries. In this study we calculate a 
suite of stellar structure models and their non-radial pulsations to investigate why the predicted 
asymptotic g-mode £ = 1 period spacing AIIi is systematically lower than is inferred from 
Kepler field stars. We find that only models with large convective cores, such as those cal¬ 
culated with our newly proposed “maximal-overshoot” scheme, can match the average AIIi 
reported. However, we also find another possible solution that is related to the method used 
to determine AIIi: mode trapping can raise the observationally inferred AIIi well above its 
true value. Even after accounting for these two proposed resolutions to the discrepancy in av¬ 
erage AHi, models still predict more CHeB stars with low AHi ( < 270 s) than are observed. 
We establish two possible remedies for this: i) there may be a difficulty in determining AHi 
for early CHeB stars (when AHi is lowest) because of the effect that the sharp composition 
profile at the hydrogen burning shell has on the pulsations, or ii) the mass of the helium core 
at the flash is higher than predicted. Our conclusions highlight the need for the reporting of 
selection effects in asteroseismic population studies in order to safely use this information to 
constrain stellar evolution theory. 

Key words: asteroseismology — stars: evolution — stars: horizontal-branch — stars: interi¬ 
ors 


1 INTRODUCTION 

In low-mass stellar evolution the core helium burning (CHeB) 
phase is the third stage of nuclear burning - after core and shell 
hydrogen burning. In evolution calculations the core structure dur¬ 
ing this phase is highly uncertain, but it has been postulated for 
more than four decades that CHeB stars develop a zone of slow 
mixing, or “semiconvection”, beyond the fully convective core (e.g. 
Schwarzschild & Harm 1969). How this develops in models, if at 
all, depends on the treatment of convective boundaries - a major 
source of uncertainty generally in stellar models (see Section 2.3). 
Later evolution depends on the structure at the end of this phase, 
and this is typically where the results of different stellar evolution 
codes begin to diverge (see e.g. Castellani et al. 1971b; Fig. 15 in 
Paxton et al. 2013). 


* E-mail: thomas.constantino@monash.edu 


I.I Brief overview 


In this paper we make use of potent new constraints on the struc¬ 
ture of core helium burning models - mixed modes of oscillation 
detected from asteroseismology. We present calculations of non- 
radial pulsations for CHeB models evolved with a variety of mix¬ 
ing prescriptions: models with and without convective overshoot 
as well as those with a semiconvection region. We also introduce 
a new algorithm for core mixing during the CHeB phase that can 
better match the asymptotic g-mode period spacing inferred from 
asteroseismology with the use of otherwise standard physics. Fi¬ 
nally, we compute pulsation spectra at different stages of the core 
helium burning phase. In order to make this paper more accessi¬ 
ble to non-experts, we provide an extensive summary in Section 4, 
which includes clear references to previous sections that contain 
more technical detail. 
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1.2 The problem of convective boundaries 


Core helium burning stars, in their various (observational) flavours, 
may be known as subdwarf B (sdB), horizontal branch (HB), RR- 
Lyrae, red clump (RC), or secondary clump (Girardi 1999) stars. 
Most of the variation between these types of CHeB stars is due 
to differences in the mass of the convective envelope and the 
hydrogen-exhausted core beneath it (where all of the hydrogen has 
been burned to helium). The common thread between them is that 
they are ah thought to contain a central helium-burning convection 
zone that is surrounded by a helium-rich region that is not convec¬ 
tive. 

In their convective cores, CHeB stars produce carbon via the 
triple-a reaction, and oxygen via ^^C(q, 7)^®0. This burning pro¬ 
duces a growing abundance discontinuity at the formal boundary 
of the convection zone if there is no convective overshoot to induce 
mixing beyond it. This boundary is usually defined as the point 
of convective neutrality (i.e. where a displaced fluid element ex¬ 
periences no acceleration). If the Schwarzschild (1906) criterion is 
applied then this is where the radiative temperature gradient Vrad 
is equal to the adiabatic temperature gradient Vad. Therefore the 
criterion for convective stability is 


Vrad ^ Vad 5 


( 1 ) 


where 


Vad = 


/ dlnT \ 

V^lnpAd 


V rad 


f dlnr \ 

UinpAad’ 


( 2 ) 


and Vrad is the temperature gradient required for radiation to carry 
the total energy flux. In low-mass CHeB models the location of this 
boundary is unstable. The increasingly C- and 0-rich mixture in the 
convection zone is more opaque than the He-rich material just be¬ 
yond the boundary. Because of this, any mixing from convective 
overshoot (which has a sound physical basis because the boundary 
is defined only as the point of zero acceleration but where convec¬ 
tive elements still carry momentum; see Castellani et al. 1971a for 
a quantitative analysis) will increase the opacity and therefore Vrad, 
and cause the convection zone to grow. Schwarzschild & Harm 
(1969) first found that in CHeB models a partially-mixed or “semi¬ 
convection” region can then develop (see also Schwarzschild 1970; 
Paczyhski 1970). While only marginally stable according to the 
Schwarzschild criterion (Vrad ~ Vad), the semiconvective region 
is stable when the effect of the molecular weight gradient is con¬ 
sidered. This is accounted for in the Ledoux (1947) criterion for 
convective stability 

Vrad < Vad + (3) 


where 



Semiconvection regions are usually defined as stable (not convec¬ 
tive) according to the Ledoux criterion but convective (or neutral) 
according to Schwarzschild. 

In early studies, algorithms were developed that produce what 
we shall refer to as “classical semiconvection”, where the composi¬ 
tion is adjusted to produce Vrad/Vad = 1, which results in a smooth 
abundance profile (e.g., Simpson 1971; Robertson & Faulkner 
1972; Faulkner & Cannon 1973). Even if there is no explicit pro¬ 
cess for allowing semiconvection, a similar chemical profile and 
temperature gradient is produced from localized mixing episodes in 
evolution sequences that have instantaneous mixing in convection 


zones and an overshooting prescription that allows mixing beyond 
the Schwarzschild boundary (Lattanzio 1986; Caloi & Mazzitelli 
1990). The most obvious difference is that the local mixing events 
leave behind numerous small composition discontinuities. By the 
end of core helium burning, both of these schemes permit the par¬ 
tially mixed region to grow to such an extent that its enclosed mass 
is around double that of the convective core (Figure 2b). Typically, 
models initially experience a period of rapid growth of the con¬ 
vective core. The expansion rate of the convective core then slows 
as a result of the emergence and subsequent growth of a partially 
mixed region. The semiconvection region in the 1.5 Mq model 
from Faulkner & Cannon (1973), for example, appears when the 
central helium mass fraction has reduced to 0.75 (from 0.978 ini¬ 
tially), at which time the convective core growth rate is halved. 

The total mass of helium that bums during the CHeB phase 
differs from code to code. The principal reason is mixing: specif¬ 
ically whether the criterion for convection is Schwarzschild or 
Ledoux, and whether convective overshoot or a scheme for semi¬ 
convection is applied. Recently, Gabriel et al. (2014) also high¬ 
lighted the significance of whether the location of a convective 
boundary is determined from inside or outside the convection zone. 
The greater opacity of the products of helium burning means that in 
this phase, numerical subtleties such as these have a compounding 
effect on the evolution. 


1.3 Classical constraints from globular clusters 

Despite making a vast difference to the evolution of the interior, the 
core mixing has little immediate effect on the conditions at the sur¬ 
face. By controlling the amount of helium that is burnt, the mixing 
scheme does, however, affect the CHeB and early-asymptotic giant 
branch (AGB) lifetimes. Empirically, the lifetime of various phases 
of evolution can be inferred from star counts in globular clusters. 
This is because they have large and (relatively) homogeneous stel¬ 
lar populations. By using the so-called R-method on a sample of 
15 globular clusters (i.e. determining Ri, the ratio of AGB to red 
giant branch stars), Buzzoni et al. (1983) found indications for the 
existence of a fully developed semiconvective zone. 

Late in the CHeB phase, models can also show the phe¬ 
nomenon of “core breathing pulses”, the rapid growth in the mass 
of the convective core when the central helium abundance is 
very low (Sweigart & Demarque 1973; Castellani et al. 1985). De¬ 
spite their emergence in stellar evolution calculations, Caputo et al. 
(1989) and Cassisi et al. (2001) contend that evidence from star 
counts in globular clusters discredits the existence of core breathing 
pulses, because they further prolong the HB lifetime and shorten 
the early-AGB. This conflict between theoretical predictions and 
observations exposes the uncertainty of stellar models during the 
CHeB phase. 


1.4 Asteroseismology of CHeB stars 

The study of asteroseismology promises a unique chance to con¬ 
strain CHeB models. The long time series observations from the 
CoRoT and Kepler missions have yielded unprecedented poten¬ 
tial for red-giant asteroseismology. Solar-like oscillations have 
now been detected in more than 13,000 giants in the Kepler field 
(Stello et al. 2013). 

Bedding et al. (2010) first detected mixed modes in the surface 
oscillations of red giants in the Kepler field. These propagate as 
acoustic modes in the convective envelope and as gravity modes in 
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the radiative core (Aizenman et al. 1977). Crucially, the observed 
period spacing AP of the dipole (£ = 1) modes is thought to 
provide a lower bound on the asymptotic ‘pure g-mode’ spacing 
AIIi (Beck et al. 2011). Bedding et al. (2011) showed that this pe¬ 
riod spacing can reliably distinguish CHeB stars from photometri¬ 
cally similar, but shell-hydrogen burning, red giant branch (RGB) 
stars. This is possible because the mixed-mode period spacing is 
sensitive to the conditions in the core, which change substantially 
between the RGB and CHeB phases. More recently, Mosser et al. 
(2012) have developed a method to infer AHi from the relatively 
small fraction of mixed modes that are detectable. The CHeB stars 
for which they reported AHi mostly have asteroseismic scaling- 
relation-determined masses of 0.8 < M/Mq <2.6, while a hand¬ 
ful have masses up to M = 3.4 M©. The metallicity ([M/H]) dis¬ 
tribution of the stars in the latest (and larger) core helium burn¬ 
ing sample with AHi determinations from Mosser et al. (2014) is 
strongly peaked around the solar value (determined from the stars 
also in the APOKASC catalogue; Pinsonneault et al. 2014). 

Measurement of AHi is a particularly useful diagnostic be¬ 
cause it depends only on the Brunt-Vaisala frequency N, which is 
easily computed from a given stellar structure. Specifically, N is 
the frequency of oscillation that an adiabatically displaced mass el¬ 
ement will undergo due to buoyancy forces. In a convective region, 
displaced elements are buoyantly unstable, therefore such oscilla¬ 
tions cannot occur, and gravity waves are damped. In the asymp¬ 
totic limit, the gravity mode period spacing is 


AUe 




-1 



(5) 


where the integral is over the region with > 0 , £ is the spherical 
harmonic degree. 


and g is the local gravitational acceleration (Dziembowski 1977). If 
> 0 then A is real, which is equivalent to the Ledoux criterion 
for convective stability 


where 


d Inp 1 
d In p ^ Ti ’ 


Ti 


/ 91np\ 


(7) 

(8) 


and the subscript ad signifies an adiabatic change. Another key ob¬ 
servable from asteroseismology is the large frequency separation, 
whose asymptotic value Av can be computed as 



where c is the adiabatic sound speed (Vandakurov 1967; Tassoul 
1980). Under a homologous transformation this scales with the 
square root of the mean stellar density, which is a finding that is 
supported by models (Ulrich 1986). 

If Av and AP (or AHi) are both determined for a star, it can 
be placed on the Av—AP (or A;z —AHi) diagram. When this pro¬ 
cedure is performed for the Kepler field stars, two distinct groups 
are found, comprising the CHeB and RGB stars respectively (e.g.. 
Fig. 3 in Mosser et al. 2012 and Fig. 1 in Mosser et al. 2014). The 
most striking feature of the Ai^ — AHi diagram is how tightly most 
of the low-mass (M/M 0 < 1 . 75 ) CHeB stars are grouped, with 


250 < AHi (s) < 340 and Av ~ 4 MHz. To date, however, CHeB 
models have been unable to properly match the AHi inferred from 
the observations (e.g. Figure 1). 

Montalban et al. (2013) identified a linear dependence of AHi 
on the radius of the convective core for CHeB models. They also 
noticed how the AHi dependence on stellar mass is similar to that 
for Mhc, and emphasized a linear dependence of AP between “ob¬ 
servable” modes on Mhc for models massive enough to avoid the 
degenerate ignition of helium (M> 2 . 2 M 0 ). Additionally, they 
suggested that a model with a semiconvection zone will have a 
lower AHi than a model with an identically sized convective core 
but without semiconvection. Their comparison between AHi from 
models (of around solar composition and 0.7 — 3.0 M 0 computed 
with the ATON evolution code; Ventura et al. 2008) and the observa¬ 
tions reported by Mosser et al. (2012) reveals a general offset, with 
the theoretical AHi lower than observed (their Fig. 7). This offset 
is also evident from models computed with MESA (Bildsten et al. 
2012; Stelloetal. 2013) and the Monash University code MON¬ 
STAR (Constantino et al. 2014a), making it apparent in at least three 
independent evolution codes. It is not as obvious that this offset 
exists for more massive, higher-Ai^ models. However, the higher- 
mass models without overshoot by Montalban et al. (2013) do not 
match the whole observed spread of AHi. Those models have 
roughly 160 < AHi (s) < 230 compared to 145 < AHi (s) < 300 
observed. Interestingly, for the 1.5 M 0 model at least, it appears 
that convective overshoot during CHeB considerably increases 
AHi (by around 50 s). 

Recently it has been shown that additional diagnostic infor¬ 
mation about mixing events may be obtained from the effect that 
resulting sharp features in the buoyancy frequency have on the ob¬ 
served mode frequencies (Cunha et al. 2015). Such features can 
arise from composition discontinuities left by first dredge-up dur¬ 
ing the RGB evolution or, as we discuss in this paper, from those 
that may arise in the CHeB phase. 


1.5 The core mass at the flash 

Sweigart & Gross (1976) showed that the zero-age horizontal 
branch (ZAHB) convective core mass depends predominantly on 
the H-exhausted core mass Mue (and is insensitive to composi¬ 
tion and total mass). Since there is a close relationship between 
radius of the convective core (and therefore its mass) and AHi 
(Montalban et al. 2013), there must also be one between Mhc and 
AHi. The mass of the H-exhausted core is therefore a crucial quan¬ 
tity for the seismology of CHeB stars. 

In low-mass CHeB models (M< 2 . 2 M 0 ), helium ignition 
occurs under degenerate conditions, triggering a thermal runaway 
known as the core flash. The minimum Mhc required for this 
ignition is fairly constant over a wide range of stellar masses 
and only decreases slightly with increasing metallicity or helium. 
Salaris & Cassisi (2005) show, for instance, that there is a 0.03 M 0 
difference in Mhc between Z = 0.02, Y = 0.273 and Z = 10“^, 
Y = 0.245 solar-mass models. Catelan et al. (1996) explored 
the uncertainties in the physics in stellar models that could influ¬ 
ence Mue, including neutrino losses, rotation, conductive opacity, 
coulomb effects on the equation of state, reaction rates and screen¬ 
ing factors, and element diffusion. They constrained the possible 
core mass increase compared with standard models to AMhc = 
0.01 ±0.015 M 0 . The best observational constraint on Mne comes 
from comparisons with globular clusters, but additional factors 
must be considered there as well, such as the composition, MLT 
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mixing length, bolometric corrections, and distance determination 
(Viaux etal. 2013). 


1.6 Clues from subdwarf B stars 

In contrast to the more massive red clump stars, the core mass of 
sdB stars is less uncertain because of the very thin hydrogen en¬ 
velope. The mass distribution of sdB stars, determined from as- 
teroseismology and eclipsing binaries, is peaked at ~ 0.47 Mq 
(V an Grootel et al. 2013a). This corresponds closely to the canoni¬ 
cal core mass at the core flash. What we may be able to learn from 
this mass distribution though is dependent on understanding the 
conditions under which helium ignition occurs. When a sample of 
sdB stars is plotted in the log g — Tcs plane there is evidence for 
two distinct groups (Green et al. 2008), perhaps suggesting differ¬ 
ent formation channels (Van Grootel et al. 2013a). Proposed binary 
mechanisms include common-envelope ejection, stable Roche lobe 
overflow and double helium white dwarf mergers (Han et al. 2002). 
If an sdB star is formed via the stable Roche lobe overflow channel 
then it is likely that its mass is close to the minimum H-exhausted 
core mass required for helium ignition. Even if we cannot reliably 
deduce the H-exhausted core mass of other CHeB stars from the 
empirical mass of sdB stars, their observed pulsation properties can 
still serve as a useful constraint on the physics during the CHeB 
phase. 

Space based observations of g-mode sdB pulsators have 
proven to be superior to earlier efforts from the ground 
(Charpinet et al. 2010). Structural properties, such as core and 
envelope mass and central helium abundance, of a handful of 
such stars have now been estimated (Van Grootel et al. 2010a,c,b; 
Charpinet et al. 2011). In their method they determine these quan¬ 
tities by finding a least-squares fit between the observed and the 
theoretical periods from models with different structural parame¬ 
ters (Brassard & Fontaine 2008). It has been reported, for example, 
that the sdB star KPD 1943-1-4058 has a larger mixed core (defined 
as the mass interior to the unmixed He-rich radiative zone) than 
models by Dorman et al. (1993) which do not include convective 
overshooting. It is unclear, however, whether the determined struc¬ 
ture includes a semiconvection region. 

Reed et al. (2011) found a period spacing range of 231 < 
AHi (s) < 271 with an average of 254 s for 13 g-mode pulsating 
sdB stars observed by Kepler and another by CoRoT. This period 
spacing is clearly lower than the more massive CHeB stars in the 
Mosser et al. (2012) sample, which have reported AHi typically 
around 300 s. 


1.7 The challenges and potential of CHeB asteroseismology 

The pulsations in CHeB models can be far more complex than they 
are in RGB models. Importantly for the propagation of g-modes, 
there is only a single radiative zone in RGB stars, whose structure 
is well understood, and apart from the discontinuity left by first 
dredge-up, it contains only smooth variations in chemical compo¬ 
sition. In contrast, the structure of CHeB models is sensitive to the 
treatment of convective boundaries in the core as well as to prior 
evolution, especially the core flash. They may contain multiple con¬ 
vection zones and composition discontinuities. The deficit in our 
understanding of the mixing during core helium burning spans a 
broad mass range, from M>0.47Mq to at least M ~ 20 Mq 
(Danger 1991). Part of the reason for this uncertainty is that until 
now the core mixing has been hidden from view. 



MJMo 

Figure 1. Compaiison between AIIi infeired from observations of likely 
CHeB stars with seismic mass determinations (black circles; Mosser et al. 
2014) and the average computed from CHeB models with standard over¬ 
shoot (cyan line; with markers showing each calculation). The line of best 
fit for observations (red dashes) follows the mode of the AHi distribution. 
The shaded area gives the range of AHi in which the models spend 95 per 
cent of their CHeB lifetime. 

The detailed study of pulsations in models of CHeB stars is 
imperative if we are to properly interpret asteroseismic observa¬ 
tions and gain understanding about the behaviour of convection 
near the boundary of the convective core. In this paper we ad¬ 
dress this need by analysing the non-radial pulsations in a range 
of CHeB models with disparate internal structures stemming from 
different treatments of convective boundaries. Any insights about 
this mixing will also have implications for the treatment of convec¬ 
tive boundaries in stellar models more generally. 

2 METHODS 

2.1 Evolution code 

In this study our evolution models are computed with MON¬ 
STAR, the Monash University stellar structure code, which 
has been described in detail previously (e.g., Lattanzio 1986; 
Campbell & Lattanzio 2008; Constantino et al. 2014b). Unless 
stated otherwise, our models are IMq, solar metallicity 
(Asplund et al. 2009), with initial helium mass fraction Y — 0.278. 

2.2 Pulsation analysis 

The models we use for pulsation calculations are usually mid-way 
through CHeB and have central helium abundance Y = 0.4 or 
Y = 0.5. The pulsations are computed with the Aarhus adiabatic 
oscillation package ADIPLS (Christensen-Dalsgaard 2008). In this 
study we restrict our analysis of non-radial modes to the t = 1 case. 
Each structure model we present is converged and in hydrostatic 
equilibrium, and the Brunt-Vaisala frequency is calculated directly 
from P, p, r, and Ei according to Equation 6. 
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Figure 3. Internal properties of the 1 Mq models with four different mix¬ 
ing prescriptions when Y = 0.4 in the centre. The helium mass fraction 
Y, ratio of temperature gradients Vj-ad/Vad. and BiTJUt-Vaisala frequency 
are shown. The four mixing prescriptions are no overshoot (Section 2.3.1; 
black), standard overshoot (Section 2.3.2; cyan), classical semiconvection 
(Section 2.3.3; orange), and maximal overshoot (Section 2.3.4; magenta). 
This colour scheme is used for mixing comparisons throughout this paper. 


Figure 2. Evolution of internal helium abundance during CHeB with four 
different mixing prescriptions. Each panel shows the profile at five dif¬ 
ferent times. The mixing prescriptions are, from top to bottom, no over¬ 
shoot (Section 2.3.1; black), standard (pure Schwarzschild) overshoot (Sec¬ 
tion 2.3.2; cyan), semiconvection (Section 2.3.3; orange), and maximal 
overshoot (Section 2.3.4; magenta). 


The estimates of the frequency of maximum power i/max for 
the models in this paper use the assumption that it scales with the 
acoustic cut-off frequency (Brown et al. 1991), and therefore that 


which is defined as 


Y; [e+m+mi]pr^dr 


( 11 ) 


where i?s is the radius at the outermost point, ri is the location of 
the innermost mesh point, and and are the radial and horizon¬ 
tal displacement eigenfunctions, respectively, which are both func¬ 
tions of r. This is a measure of kinetic energy of a mode relative to 
the radial displacement at the surface. In the plots of eigenfunctions 
we show the scaled horizontal displacement 


( 10 ) 

where g is the surface gravity, qq = 2.74 x 10“* cm s“^, r/max,© = 
3.1 mHz, and Teff,© = 5778 K. 

We obtain frequencies of individual modes from the pulsa¬ 
tion calculations and hence the actual period spacing AP between 
modes of adjacent order. We present these results by showing AP 
as a function of frequency. In several examples we use the period 
echelle diagram, where the mode frequency is plotted against the 
mode period, modulo some period spacing APech- This is used be¬ 
cause g-modes tend to be approximately equally spaced in period 
(the asymptotic limit is given in Equation 5). This is important be¬ 
cause in the method developed by Mosser et al. (2012), AHi corre¬ 
sponds to the APcch which produces a regular pattern in the echelle 
diagram. Thus it allows us to predict the value of AHi that would 
be inferred from observations of our theoretical models. 

We report the mode inertia E from the pulsation calculations 


g_ 

90 


Teff,© 


Tef 


i^max,© ? 


y2 


£{£+!) 

R 


Ch, 


( 12 ) 


where R is the photospheric radius and is scaled so that ^r/R ~ 
1 at the surface. 


2.3 Core mixing schemes 

2.3.1 Models without convective overshoot 

In the stellar models without convective overshoot the 
Schwarzschild criterion for convection is strictly applied. 
The resulting internal structure is shown in Figure 2a and Figure 3. 
In this case a convection zone may only grow (in mass) if the 
conditions change in a radiative zone so that V,ad > Vad- In the 
CHeB phase this will give the smallest possible convective core 
(at least for models with the Schwarzschild criterion) because 
the region outside the core is close to convective neutrality (see 
Figure 3). In this study we do not compute any models using the 
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Ledoux criterion for convection, but we note that we do not expect 
it to make any difference to our models either with or without 
convective overshoot, if properly implemented. This is because 
in our models without overshoot the convective boundary hardly 
moves (by less than 0.001 M© in the 1 Mq solar-metallicity run; 
Figure 2), so any restriction in growth due to the composition 
gradient would be insignificant. In the overshooting models, 
the mixing beyond the boundary tends to erase any stabilizing 
composition gradients and therefore reduces the Ledoux criterion 
to the Schwarzschild one. 


2.3.2 Models with convective overshoot 

We have implemented overshooting (Campbell & Lattanzio 2008) 
with an exponential decay in the diffusion coefficient according to 
the scheme proposed by Herwig et al. (1997). This is expressed as 

Dos{z) = Doe "v* , ( 13 ) 

where Dos(z) is the diffusion coefficient at distance z from the 
convective boundary and Do is the diffusion coefficient just inside 
the boundary. H, is the “velocity scale height” defined as 

TTv = /osTTp, (14) 

where is the pressure scale height, and we have chosen /os = 
0.001. We refer to this as “standard overshoot”, but the exact value 
of /os is not important because our models are insensitive to the 
formulation of convective overshoot. The resulting internal struc¬ 
ture is shown in Figure 2b and Figure 3. 

Our models with convective overshoot and time-dependent 
mixing of chemical species (Campbell & Lattanzio 2008) evolve 
similarly to those with instant mixing and the search for convec¬ 
tive neutrality (Lattanzio 1986). In the latter method, the convective 
boundary is found by testing whether mixing at the Schwarzschild 
boundary would cause the next radiative zone to become convec¬ 
tive, while the former always mixes beyond the Schwarzschild 
boundary (i.e. without a test). The outcome is similar because at 
these conditions C and O are so much more opaque than He, so 
this mixing usually results in the radiative zone adjacent to the 
Schwarzchild boundary at the outer edge of the convective core be¬ 
coming unstable to convection. This feedback contrasts with other 
phases in evolution when the extent of mixing is dependent on the 
overshooting distance because the resulting mixing does not alter 
the location of the Schwarzschild boundary. 

In CHeB models, any overshoot tends to grow the core enough 
so that a minimum in Vrad/Vad appears in the convection zone; see 
the magenta line in Figure 3b for an example. Eventually this min¬ 
imum falls below unity and the convection zone splits. This pro¬ 
cess continuously repeats, leaving behind the characteristic stepped 
abundance profile seen in Figure 2b. If properly resolved, the par¬ 
tially mixed region created beyond the convective core by over¬ 
shooting will have temperature gradient Vrad/Vad ~ 1, resembling 
semiconvection (see Section 2.3.3). 


2.3.3 Models with semiconvection 

We have developed a simple method to mimic the structure that is 
found using semiconvection routines. We do this by allowing slow 
mixing in sub-adiabatic conditions. Specifically, we set a mixing 
rate that depends only on how close a zone is to being convective 
according to the Schwarzschild criterion (which neglects the stabi¬ 
lizing effect of any composition gradients). If Vrad < Vad then we 


set the diffusion coefficient D according to 

log D = log D' - Cl (1 - Vrad/Vad), (15) 


but also specify a maximum gradient so that 


dlnZ) 

dlnp 


< C2, 


( 16 ) 


where D', ci, and C 2 are constants that are chosen at discretion 
and varied with experience (cf. Iben 1967). In this study we use 
log D' — 10, Cl = 100, and C 2 = 90 to give D in units of cm^s“^. 
The structure that is produced by this scheme is shown in Figure 2c. 
This differs slightly from previous routines that have a zone with 
exactly Vrad = Vad and a smooth composition profile which ends 
with a discontinuity (herein the “classical” semiconvection struc¬ 
ture). An example of the classical structure is shown in Figure 3. 
The most obvious difference produced by our routine is that the 
composition is everywhere smoothly varying, i.e. there is no dis¬ 
continuity at the outer boundary of the partially mixed region. 

We also use a different method to construct classical semicon¬ 
vection models. For a given central composition we artificially in¬ 
crease the mass fraction of the convective core until there is a min¬ 
imum in Vrad/Vad inside the convection zone with Vrad/Vad = 1 
(the structure during this intermediate step is identical to that of 
the maximal-overshoot models shown in Figure 3 and discussed 
below). The location of this minimum is the first guess for the 
boundary between the convection and semiconvection zones, while 
the edge of the convective region becomes the outer boundary of 
the semiconvection zone. We then adjust the composition between 
these two points until this region has Vrad/Vad = 1 everywhere 
(and also make small adjustments to the location of the boundaries 
if needed). This contrasts with the method of Robertson & Faulkner 
(1972) where, during each time step, the composition changes due 
to nuclear burning and then mixing proceeds outwards from the 
centre, point by point, to give exactly Vrad = Vad- 


2.3.4 Models with maximal overshoot 

If the helium-burning convective core is large enough, it will con¬ 
tain within it a minimum Vrad/Vad, such as that shown by the ma¬ 
genta line in Figure 3b. Further core growth will continue to reduce 
the value of this minimum until it reaches Vrad/Vad < 1, which 
splits the convection zone into two. This is avoided in our newly de¬ 
veloped “maximal-overshoot” scheme by making convective over¬ 
shoot dependent on the value of this minimum, so that the core 
growth slows (and then can stop) if the convective core is close to 
splitting. In its present ad hoc implementation the amount of mass 
beyond the convective boundary that is mixed each time step is pro¬ 
portional to the minimum of Vrad/Vad — 1 in the convection zone, 
and overshoot is stopped if that minimum falls below 0.002. This 
ensures that the model attains the largest possible convective core 
throughout the evolution. This structure is shown in Figure 2d and 
Figure 3. We do not propose a physical justification for achieving 
this exact structure. Instead we use it as a comparison to standard 
models which is interesting because of its extreme core size and 
the effect on AHi. Finally, we note that although our maximal- 
overshoot models are generated by a different mechanism, their 
structure is similar to some earlier models with large fully mixed 
cores (e.g., Bressan et al. 1986; Straniero et al. 2003). 
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Figure 4. Left panel: Kippenhahn plot of the core flash phase of a 1 Mq 
solar-metallicity model. Right panel: schematic Kippenhahn plot of sub¬ 
sequent quiescent CHeB evolution. Shaded areas (grey) denote convection 
zones. The upper curve (blue) and lower curve (magenta) denote the re¬ 
spective positions of maximum H- and He-buming luminosity. In the right 
panel the grey stripes show the region where partial mixing or semiconvec¬ 
tion may occur. 



Afcc/Mo Rcc/Rq Mke/Mo 


Figure 5. Dependence of the asymptotic g-mode period spacing AIIi on 
core properties (mass of the convective core Mcc, radius of the convec¬ 
tive core i?cc, and mass of the H-exhausted core Mhc) for models with the 
smallest (black) and largest (magenta) possible convective core, artificially 
constructed by varying Mue. The models are IMq, solar-composition, 
and consist of a chemically homogeneous convection zone beneath the H- 
burning shell. Each model has central helium Y = 0.5. 


2.4 Composition smoothing 

The period spacing pattern in CHeB models taken directly from the 
evolution code tends to be inconsistent with observations (e.g. Fig. 

2 in Mosser et al. 2012). In Section 3.3 we analyse this in detail and 
show that it is primarily a relic of the burning during the core flash 
phase. Therefore, in a number of our ZAHB models we remove the 
composition discontinuities between the H-burning shell and the 
convective core that this burning produces. We have good reason 
to do this; the chemical profile that the core flash leaves behind 
is dependent on unknown factors such as the extent of convective 
overshoot and mixing at boundaries, the extent of burning during 
each episode of convection (or “subfiash”) as burning progresses 
inward, and the number of these subflashes. We remove this fea¬ 
ture either by artificially resetting the composition to the mixture 
that existed immediately prior to the core flash or by smoothing the 
composition over a larger interval in mass. In the latter method we 
set the mass fraction Xi of species i according to 

Xi 1 -f Xi 2 Xi 2 — Xi 1 
'2 ’ + ’2 ’ 

where Am is the mass over which the composition is smoothed, 
centred at mq, and and Ai ,2 denote the interior and exte¬ 
rior compositions. After modifying the composition using either 
method we then evolve the model in the evolution code to allow 
it to return to hydrostatic equilibrium before computing the pulsa¬ 
tions. 

3 RESULTS 

3.1 Overview of the CHeB structure 

We show a schematic of the evolution of the internal structure of 
a CHeB model in Figure 4. The profile of the Brunt-Vaisala fre¬ 
quency that develops is crucially dependent on the mixing scheme 
used at the boundary of the convective core, which is evident from 
the difference between the models in Figure 3c. The significant fea¬ 
tures affecting the Brunt-Vaisala frequency (proceeding outward 
from the centre) are as follows. 


(i) A fully mixed convective core that is Ledoux and 
Schwarzschild unstable {N ~ 0). 

(ii) A region that may surround the convective core in which 
material is slowly mixed (see e.g. the grey stripes in Fig¬ 
ure 4). Depending on the mixing scheme there can emerge 
zones with a stabilizing chemical gradient (Schwarzschild 
marginally stable, Ledoux stable, and A 0) or regions 
which are convective and well-mixed (N ~ 0). The erratic 
nature of overshoot can create an irregular N profile in this 
region that constantly evolves. In this study, such a region 
only emerges in the standard-overshoot (Section 2.3.2) and 
semiconvection (Section 2.3.3) models. 

(iii) A helium-rich radiative region with A > 0. In models in 
which helium ignition begins with the core flash there will be 
composition gradients between the (fully or partially mixed) 
core and the H-burning shell (blue line in Figure 4). Only 
a small fraction of the helium bums in the core flash phase 
(around 3 per cent in our 1 Mq models) but this is enough 
to cause detectable spikes in A from the molecular weight 
gradients formed at the boundaries of flash and subflash con¬ 
vection zones, e.g. near r = 4.2 x 10® cm in Figure 9. In our 
models the largest spike is caused by the burning in the initial 
core flash (closest to the H-buming shell). 

(iv) The H-burning shell, which is strongly stable due to the 
molecular weight gradient (A ^ 0; blue line in Figure 4). 

(v) A radiative zone below the convective envelope (or the sur¬ 
face if the star is not massive enough to have a convective 
envelope). A decreases monotonically until the convective 
envelope, where A ~ 0. The convective envelope extends 
all the way to the surface (or close enough for the pulsation 
calculations we perform in this study). 

The radiative region between the convective envelope and the con¬ 
vective core is where g-modes propagate. Importantly, this includes 
any partially mixed or semiconvection region surrounding the con¬ 
vective core (e.g., the grey stripes in Figure 4). 


m — mo 


Am 


(17) 
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Figure 6. AIIi dependence on the mass, radius and minimum of 
Vrad/Vrad of the fully mixed core Mcore- These models are constructed 
by varying Mcore while keeping Mfje = 0.488 Mq and the core convec¬ 
tion zone composition constant with Y = 0.5. Each model is 1 Mq and 
solar metallicity. 


3.2 AIIi dependence on bulk core properties 

We have tested the dependence of AIIi on the fundamental proper¬ 
ties of the core - the mass of the H-exhausted core Mhc (Figure 5) 
and the mass of the fully mixed core Mcore (Figure 6). 

In Figure 5 we increase Mhc through H-burning then con¬ 
struct models with the smallest (black markers) and largest (ma¬ 
genta markers) possible convective cores for a fixed central compo¬ 
sition of y = 0.5. In this test the mass of the convective core Mcc 
(Figure 5a) and radius of the convective core Rcc (Figure 5b) are 
dependent variables. The smallest convective core is found by ex¬ 
tending it only until the region adjacent to the convection zone (i.e. 
the radiative side) is marginally stable to convection, which gives 
the same structure as the “no overshoot” models in this study. In 
contrast, the models with the largest core are constructed by extend¬ 
ing the convective core as far as possible so that the entire region 
within it remains convectively unstable, which is the same as for 
the maximal-overshoot models (Figure 3, Section 3.4.4). We find 
that in both the smallest and largest convective core cases Mcc, i?cc, 
and AIIi are linearly dependent on Mhc. Montalban et al. (2013) 
have already highlighted the linear dependence of AIIi on i?cc in 
low-mass CHeB models, as well as the importance of Mhc, includ¬ 
ing a linear relationship between AIIi and Mhc for more massive 
models. The difference in AIIi between the smallest and largest 
core cases is 45 s. It is evident from the fact that two values of i?cc 
can correspond to the same value of AIIi that there is not a single 
linear dependence of AIIi on i?cc (Figure 5b). In this case the rela¬ 
tionship also depends on how the convective boundary is defined. 

In Figure 6 we show the effect of artificially changing the 
mass of the homogeneous (fully mixed) region in the core, Mcore, 
while keeping the central composition and the H-exhausted core 
mass constant. Note that Mcore differs from Mcc in that there is 
no requirement that the entire region enclosed by Mcore is convec¬ 
tive according to the Schwarzschild criterion. The peak in AHi 
occurs exactly when the convective zone is as large as possible 
(Mcore = 0.22 Mq). This clearly demonstrates that further extend¬ 
ing the fully mixed core (e.g. as a result of overshoot; such as the 
“high overshoot” model from Straniero et al. 2003) does not con¬ 
tinue to increase AHi when part of it becomes stable to convection 
according to Schwarzschild, allowing g-modes to propagate. This 
is of interest because it is not unreasonable to imagine that convec¬ 
tive overshoot could allow the composition of two nearby convec¬ 


tion zones to remain homogeneous. In this example, models with 
0.22 < Mcorc/MQ < 0.33 have two separate convection zones. 

In Figure 7 we compare the £ = 1 period spacing for mod¬ 
els with fully mixed cores of different sizes: one with a fully 
mixed core mass of 0.215 Mq (black dashes) and another other¬ 
wise identical model, with a fully mixed core mass of 0.255 Mq 
(orange dashes), which is too large to be convective throughout 
and therefore gives rise to an additional radiative region (1.7 x 
10®cm<r<2.3 x 10^ cm in Figure 7). The asymptotic period 
spacing (from integrating over the entire structure according to 
Equation 5) is slightly lower for the model with the larger core 
(by 3 s). Some modes in this model are very closely spaced in pe¬ 
riod because they are trapped in the additional radiative region (see 
also Sections 3.4.2 and 3.4.3), whereas most pairs of modes have 
AP > AHi. If the mode periods are plotted in the echelle dia¬ 
gram (see Section 2.2), the period spacing required for a regular 
pattern, APech, is 19 s higher than the asymptotic value. The reason 
for this can be understood by considering the local buoyancy radius 
n“^(r) described by Miglio et al. (2008). At a given point r, this 
is defined by 

n’^W= f ^dr', (18) 

J rg 

where ro is the radius at the edge of the convective core (r ~ 1.7 x 
10® cm in Figure 7). This gives the contribution to the integral in 
Equation 5 from the region enclosed by the point at radius r. The 
total buoyancy radius is the same integral evaluated over the entire 
g-mode propagation zone. 

In the model with the larger fully mixed core (orange dashes 
in Eigure 7) the interior buoyancy cavity accounts for 7.8 per cent 
of the total buoyancy radius, which corresponds to the difference 
between AHi and APech (which are 284 s compared to 307 s re¬ 
spectively). If we exclude this interior cavity from the calculation 
of AHi we get almost exactly AHi = APech (308 s compared to 
307 s). This divergence between AHi and APech is also apparent 
for the model with a sharp composition profile in Eigure 7 (in blue), 
demonstrating that it does not depend on the composition profile at 
the edge of the fully mixed core. We therefore expect that the exis¬ 
tence of a second radiative zone would generally cause the observa- 
tionally inferred value of AHi (using the method of Mosser et al. 
2012) to increase above its theoretical value (computed by integrat¬ 
ing over the entire structure). This example highlights the possibil¬ 
ity that a difficulty in accurately determining AHi from observa¬ 
tions may contribute to its apparent discrepancy with predictions 
from models. This phenomenon is discussed in more detail for a 
related example in Section 3.4.2. 

Eigure 8 shows the evolution of the internal structure of mod¬ 
els with the four different mixing schemes. The size (in mass and 
radius) of the convective core in the semiconvection and standard 
overshoot sequences is similar throughout the evolution, except 
when overshooting permits core breathing pulses near core helium 
exhaustion. This explains the similarity in AHi evolution. In both 
the semiconvection and standard-overshoot sequences almost all of 
the growth in the mass of the convective core occurs during the first 
20Myr (Eigure 8c). Subsequently, helium is transported into the 
core by the expansion of the partially mixed region. Interestingly, 
the rate of depletion of helium in the core is exactly the same for 
the maximal-overshoot and standard-overshoot runs until the final 
(and largest) core breathing pulse extends the standard-overshoot 
model CHeB lifetime. 

At the beginning of CHeB, each of the four standard se- 
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quences (solid lines) shows a decrease in AIIi (Figure 8a). This 
can be attributed to the softening of the steep composition gradient 
at the H-burning shell. This is further discussed in Section 3.3 and 
its effect on AIIi is also explained for an analogous case in Sec¬ 
tion 3.5. After hydrogen burning resumes over the entire shell, the 
evolution of AIIi closely tracks the radius of the convective core, 
which has been shown by Montalban et al. (2013). 

In dashed lines in Figure 8a we also show additional sequences 
with standard and maximal overshoot that result from enlarging 
Mue by AMue = 0.025 Mq at the beginning of CHeB. This was 
achieved by delaying helium ignition through an ad hoc increase 
to the neutrino emission rate during the RGB phase. This increases 
the average AIIi during CHeB by 18 s for the standard-overshoot 
sequence and 11s for the maximal-overshoot case. Most signifi¬ 
cantly, it increases AHi by around 20 s early in the CHeB phase for 
both sequences. The faster rate of helium burning resulting from the 
larger Mhh exhausts the fuel earlier, shortening the CHeB lifetime 
by around 25 per cent for both mixing schemes. 

As core helium burning progresses, the convective core be¬ 
comes increasingly C- and 0-rich (Figure 8d) and consequently 
more dense. This causes the convective core radius to decrease 
(Figure 8b), even when its mass does not and irrespective of the 
mixing scheme (Figure 8c). It is evident from the decrease in AHi 
towards the end of CHeB for every sequence shown in Figure 8a 
that AHi is more closely dependent on convective core radius than 
mass. The final composition of the degenerate C-O core is another 
potential diagnostic for mixing, but the range covered by these four 
different schemes is small (oxygen varies by around 15 per cent) 
and the situation is further complicated by reaction rate uncertain¬ 
ties (Straniero et al. 2003). 

In this section we have explored how differences in the stel¬ 
lar structure affect AHi. Even among the models without a semi¬ 
convection or partially mixed zone, AHi depends on a number of 
factors: Rcc, Mhc, and core composition. Additionally, if the chem¬ 
ically homogeneous region in the core is large enough for part of it 
to become radiative, mode trapping can cause the period spacing to 
increase above the asymptotic value. Such an effect would increase 
the AHi inferred from observations, and therefore help to explain 
why standard CHeB models do not match the average AHi for the 
Kepler field stars. 

3.3 Pulsations in early post core-flash CHeB models 

Neutrino emission from plasmon decay during the RGB phase is 
strongest at the centre, where the density is highest. This energy 
loss is enough to move the position of maximum temperature, and 
therefore He-ignition, off-centre. After helium ignition a sequence 
of several subflashes move inward until the burning reaches the 
core and quiescent CHeB begins (Figure 4). In one of the first 
studies making use of the mixed mode detection in red giants, 
Bildsten et al. (2012) found that it may be possible to identify stars 
in the core flash phase by using the fact that (between the sub¬ 
flashes) the g-mode period spacing is expected to be much lower 
than for quiescent CHeB stars, but still higher than for RGB stars. 
Their approach of studying the population in the Kepler field, if 
successful, could reveal the lifetime of the core flash phase and the 
nature of the mixing, and therefore also shed light on the stmc- 
ture before the core flash (which is dependent on neutrino losses 
for example). Here we examine the computed pulsation spectra of 
a model in the early post core-flash phase, and test the effect of the 
remaining abundance profile. 

The inward progression of convection and burning during the 




V (pHz) 



Figure 7. Upper panels: Brunt-Vaisala frequency and £ = 1 period spacing 
for three models with different sized fully mixed cores and different chem¬ 
ical profiles near the boundary. Here we have (i) a smooth boundary with a 
fully mixed core that is convective (black), and (ii) a fully mixed core too 
large to be convective with a smooth (orange) and (iii) sharp (blue) compo¬ 
sition profile at the boundary. The composition profiles are set according to 
Equation 17 with Am = 0.01 Mq (smooth profile) and Am = 10“® Mq 
( sharp profile). Lower panels: period echelle diagrams for the three models. 
Larger symbols correspond to lower mode inertia. The best fits for AP&h 
are achieved for period modulo 288 s, 307 s, and 323 s. This compares to 
the asymptotic g-mode period spacing (AIIi) of 287 s, 284 s, and 295 s, re¬ 
spectively. Each model has approximately R = 10.7Rq, Teff = 4760 K, 
and i/max = 27 /iHz. 


core flash phase has a lasting effect on the Brunt-Vaisala frequency. 
Two of these features in our early-CHeB model are visible in Fig¬ 
ure 9: the peak in A at r = 1.5 x 10® cm, due to the initial reces¬ 
sion of the convective core at the beginning of CHeB; and another at 
r = 4.2 X 10® cm, which is caused by the first episode of core-flash 
burning. However, from a seismic perspective, the dominant fea- 
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Figure 8. Evolution of the various 1 Mq models with four different treat¬ 
ments of convective boundaries during the CHeB phase. Models have the 
same colours as Figures 2 and 3. Properties shown, from top to bottom, are 
the asymptotic g-mode period spacing AIIi, the radius of the convective 
core Rcc , the mass of the H-exhausted and the convective core, and the cen¬ 
tral helium (solid line), carbon (dots) and oxygen (dashes) mass fractions. 
In the top panel dashed lines indicate models with an increased core mass 
AMhc = 0.025 Mq at the beginning of the CHeB phase. 


ture of early-CHeB models is the thinness of the H-burning shell, 
which is due to the relatively steep temperature gradient in the prior 
luminous RGB phase. In Figure 9 we confirm that the composition 
profile at the H-burning shell causes the irregular period spacing 
pattern by showing that it still exists for a model with the sharp N 
feature from the first episode of core-flash burning removed (blue 
model). 

In our solar-mass runs it takes more than 14 Myr for hydro¬ 
gen burning to smooth out the composition gradient at the inside 
of the shell at r = 4.5 x 10® cm, as shown in Figure 10, and thus 
for a regular pattern in the period spacing to emerge. In Figure 11 
we show the effect of this sharp composition gradient on the eigen¬ 
functions. After the core flash there is a sharp peak in the buoyancy 
frequency (at r = 4.3 x 10® cm in Figure 11a) which traps modes 
of consecutive radial order to very different extents (Figure 11c). 
This buoyancy peak is then slowly eroded by hydrogen burning 
(Figure 11b), and once the hydrogen burning shell completely re¬ 
activates the buoyancy peak is broad compared to the characteristic 
eigenfunction wavelength (Figure lid), leaving the period spacing 
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V (pHz) 

Figure 9. Pulsation properties of 1 M0 models with standard overshoot 
2.7 Myr after the onset of convective core helium burning taken directly 
from the evolution code (black) and with the region outside the convective 
core homogenized (blue). First panel: Brunt-Vais^a frequency N. Second 
panel: helium mass fraction Y . Third panel: 1=1 period spacing. Fourth 
panel: I = {) and £ = 1 mode inertia nonualized to the surface radial 
displacement. The model from the evolution code has Afli = 240 s, R = 
11.6 Rq, Teff = 4580 K, and Umax = 23//Hz whereas the model with 
the homogenized composition has Afli = 238 s, R = 11.5 Rq, Tgff = 
4570 K, and Umax = 24 //Hz. 

more regular. If this picture is true for real stars we anticipate diffi¬ 
culty in determining AIIi for up to 15 per cent of red clump stars 
(based on a 100 Myr CHeB lifetime) and also every star in the core 
flash phase. 

Only after hydrogen burning resumes throughout the shell 
does the effect of any composition profile left by the core flash 
become dominant (it is responsible for the difference between the 
two models in Figure 10). In Figure 12 we show that the nature 
of this chemical profile strongly affects the pulsations. We com- 
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pare two chemical profiles that vary smoothly according to Equa¬ 
tion 17 and one that has the composition profile from the evolution 
code. It is clear in the case with the composition gradient spread 
over Am = 0.002 Mq (central panels in Figure 12) that the high- 
frequency (low radial order) modes are most sensitive to this fea¬ 
ture (blue model in Figure 12). This is because adjacent modes 
are affected differently when the characteristic eigenfunction wave¬ 
length is large enough to be comparable to the size of the peak in 
N. This causes the variation in AP between consecutive pairs of 
modes. If the composition profile is smoother, the effect of this fea¬ 
ture diminishes (e.g., in the orange model with Am = 0.01 Mq; 
left panel in Figure 12), and the period spacing resembles models 
where it is absent (e.g. the blue model in Figure 10). Conversely, 
the sharper composition profile produced by the evolution code 
(and dependent on the treatment of the convective boundary dur¬ 
ing the core flash) produces a very obvious effect on the period 
spacing. Specifically, it introduces a large mode-to-mode variation 
in AP throughout the frequency range examined (black model in 
Figure 10). This behaviour can be explained by analysing the eigen¬ 
functions. The spike in N is located at about half the buoyancy 
radius (it is at 57 per cent of the total buoyancy radius in this par¬ 
ticular case) where neighbouring eigenfunctions are separated in 
phase by about 7r/2, so it affects consecutive modes differently. In 
this case Am ~ 10“^ Mq which is small compared to the eigen¬ 
function wavelength. 

These examples serve as a note of caution when computing 
pulsations for models with steep composition gradients. The steep¬ 
ness of those resulting from the core flash phase is particularly de¬ 
pendent on convective overshoot, which is uncertain. This is the 
reason we have smoothed the composition profile after the core 
flash phase for many models shown in this paper. 

3.4 Pulsation properties for models with different mixing 
schemes 

3.4.1 Models without convective overshoot 

Our models without convective overshoot do not develop a par¬ 
tially mixed region, and experience negligible growth in the mass 
of the convective core. When the effects of core flash mixing are ex¬ 
cluded, and after H-burning broadens the shell, these models have a 
simple buoyancy profile, and consequently a simple period spacing 
pattern (Figure 13). This period spacing pattern closely resembles 
RGB models (e.g. Figure lb in Bedding el al. 2011) and observa¬ 
tions (e.g. KIC 9882316 in Fig. 1 in Mosser et al. 2012) except that 
the period spacing is higher. 

3.4.2 Models with standard overshoot 

The dominant factor in the computed pulsations of our standard- 
overshoot model shown in black in Figure 14 is the main compo¬ 
sition discontinuity left by the core flash (see Section 3.3 for the 
analysis of the effect of this discontinuity). In order to isolate the 
effect of the partially mixed region resulting from core helium burn¬ 
ing (r < 2.3 X 10® cm) we have smoothed the composition profile 
created during the core flash phase. 

In addition, we have also smoothed the chemical profile at the 
edge of the partially mixed zone (r ~ 2.3 x 10® cm), in order to 
make the period spacing slightly more regular. The resulting pe¬ 
riod spacing pattern shown in cyan in Figure 14 differs from the 
earlier models without a partially mixed region by the appearance 
of consecutive modes that are very closely spaced in period. These 
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Figure 10. Pulsation properties of 1 Mq models with standard overshoot 
14.8 Myr after the end of the core flash phase from the evolution code 
(black) and with the region outside the convective core homogenized (blue). 
Both models have AIIi = 247s. These models have R = 11.4Rq, 
Teff = 4590 K, and i/max = 24 pHz. 


have a regular dependence on radial order n and are separated by 
An « 11. These modes are also of very high inertia, and their ef¬ 
fect on AP in Figure 14 appears superimposed on the pattern pro¬ 
duced by a structure without a partially mixed region (e.g. the blue 
model in Figure 13). The reason for this is clear from Figure 15b, 
where it can be seen that these modes are “trapped” in the partially 
mixed region by the discontinuity at its boundary. The period spac¬ 
ing between other modes is affected too: the position of nodes in 
the eigenfunctions of neighbouring modes in the trapping region is 
nearly identical, so the period spacing between most of them is as 
if the interior cavity does not exist, i.e. their AP is more consistent 
with the asymptotic calculation excluding this cavity and is there¬ 
fore higher. This is demonstrated in the period echelle diagram in 
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Figure 11. Upper panels: comparison of the Brunt-Vmsala frequency N 
and helium mass fraction Y near the edge of the H-exhausted core for two 
early-CHeB models. The model represented by the solid lines is 6.8 Myr 
after the beginning of the CHeB phase and the dashed model has evolved 
for another 10 Myr. Lower panels: scaled horizontal displacement eigen¬ 
functions y 2 (defined in Equation 12) for £= I modes in c) the ear¬ 
lier, and d) the later model. The modes are of consecutive radial order 
n = —91, —92, —93 (in orange, cyan, and black, respectively), and have 
frequency of roughly 42 /rHz. These models have approximately R = 
11.3R 0 , Teff = 4570 K, and i/^ax = 24 fjHz. 

Figure 15 (see also the analogous case in Figure 7). In this exam¬ 
ple, the APcch that produces the best fit for the echelle diagram is 
315 s, which is well above the asymptotic value of 281 s. This sug¬ 
gests that inferences from asteroseismology about the theoretical 
AIIi as specified in Equation 5 may be dubious if the stars have a 
structure comparable to our standard-overshoot models. Moreover, 
the trapped modes that are responsible for this behaviour have a 
relatively small amplitude at the surface and may therefore be im¬ 
possible to detect. 

3.4.3 Models with semiconvection 

We have computed the pulsation spectra for four models with 
semiconvection-like structures. In Figure 16 we compare the classi¬ 
cal semiconvection structure to those where the abundance discon¬ 
tinuity at the outer boundary of the semiconvection zone has been 
softened. In Figure 17 we analyse the structure that is produced by 
our routine that mimics semiconvection. 

In classical semiconvection models the sharp composition gra¬ 
dient between the semiconvection zone (which is relatively C- and 
0-rich and has stabilizing composition gradient) and the He-rich 
zone produces a step in the Brunt-Vaisala frequency (at r = 
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Figure 12. Comparison of pulsation properties of the 1 Mq models with 
different composition profiles from the initial core flash. Upper panels: 
Brunt-Vaisala frequency N, helium mass fraction Y, and scaled horizontal 
displacement eigenfunctions y 2 (defined in Equation 12) for two consecu¬ 
tive i = 1 modes with u 26 /iHz. From left to right the models have 
chemical profiles artificially smoothed over 0.01 Mq (orange), 0.002 Mq 
(blue), and approximately 10“^^ Mq (from the evolution code; black). The 
former two are of the form described by Equation 17 and have approxi¬ 
mately R = 11.1 Rq, Tgff = 4750 K, and Umax = 27 //Hz while the latter 
model has R = 11.1 Rq, Tgff = 4610 K, and Umax = 25//Hz. Lower 
panels: £ = 1 mode period spacing and normalized inertia for the above 
models (same colours). 


2.2 X 10® cm in Figure 16). This causes variation in AP between 
consecutive pairs of low-frequency modes and also in the inertia of 
every second mode (Figure 16). The replacement of this discontinu¬ 
ity by a linear composition profile spread over 0.01 Mq introduces 
a second periodicity in AP of An 2 =; 7, where n is the radial or¬ 
der. This is consistent with the expression derived by Miglio et al. 
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Figure 13. Pulsation properties of 1 Mq models without convective over¬ 
shoot from the evolution code (black) and with the region outside the con¬ 
vective core homogenized (blue), and a model with maximal overshoot (ma¬ 
genta). Both models without overshoot have Afli = 252 s while the model 
with maximal overshoot has Afli = 314 s. The models have approxi¬ 
mately R = 11.1 Rq, Teff = 4610 K, and i^max = 25 //Hz. 


(2008): 

An~^, (19) 

llo 

where is the buoyancy radius at the location of the compo¬ 
sition gradient and is the total buoyancy radius (defined in 
Equation 18). This relatively smooth composition profile has only a 
small effect on period spacing and mode inertia compared with the 
case where the composition varies over just 2 x 10“^ Mq (cyan 
dashes in Figure 16). In the latter case, modes are very strongly 
trapped in the semiconvection zone (with about the same periodic¬ 
ity), which increases the period spacing between the non-trapped 



Figure 14. Pulsation properties of 1 Mq models with standard overshoot 
from the evolution code (black) and with some composition smoothing 
(cyan). The composition profile in the latter model has been smoothed near 
r = 4.2 X 10® cm by using Equation 17 with Am = 0.008 X Mq, 
which is sufficient for it to not affect the computed frequencies. The dis¬ 
continuity at the edge of the partially mixed zone (2.3 X 10® cm) has been 
smoothed with Am = 2 X 10“^ Mq. The model from the evolution code 
has AIIi = 278 s while the model with smoothing has AIIi = 281 s. 
Both models have R = 11.0 Rq, T^ff = 4600 K, and i^max = 26 pHz. 

modes to around 300 s, well above the asymptotic value of 271 s. 
This model has a regular period spacing pattern when plotted in 
the echelle diagram with APcch = 306 s, which is consistent with 
the AIIi calculation excluding the region interior to the composi¬ 
tion discontinuity. This is analogous to the model with strong mode 
trapping in Figures 14 and 15. In both cases, the low AP between 
certain pairs of modes provides the only hint that the typical AP is 
actually above the asymptotic value. One of the modes in each of 
these pairs, however, is unlikely to be detected because it has high 
inertia and is trapped in the semiconvection/partially mixed zone. 


MNRAS 000, 000-000 (0000) 
































































14 


Constantino et al. 


^1000 

£ 

100 

10 

0.1 

(N 

>1 0.0 


- 0.1 



Figure IS. Seismic properties of the standard-overshoot model in Fig¬ 
ure 14 with composition smoothing. Left panel: Brunt-vaisala frequency 
and scaled horziontal displacement eigenfunctions (where y 2 is defined in 
Equation 12) with radial order n = —130, —129, —128 (black line, or¬ 
ange line, and magenta dashes, respectively) and tz 26 /rHz near the 
abundance discontinuity at the edge of the partially mixed zone (r = 
2.3 X 10® cm in Figure 14). The eigenfunctions of the “trapped” mode 
(n = —128; magenta) has been rescaled by a factor of 0.1 for clar¬ 
ity. Right panel: Echelle diagram with APe^h = 315 s. This model has 
AIIi = 281 s, but if the calculation includes only the region exterior to the 
discontinuity then AIIi = 315 s. 


In contrast with the classical semiconvection models, the 
buoyancy spike produced by the ad hoc semiconvection scheme in 
the evolution code only weakly traps modes. This is still enough to 
clearly add a periodicity to the period spacing (with An ~ 5; Fig¬ 
ure 17). This is consistent with Equation 19 because exactly 20 per 
cent of the total buoyancy radius is contained within the partially 
mixed zone. 


3.4.4 Models with maximal overshoot 




20 30 40 50 

V (|iHz) 


The structure of the maximal-overshoot models is very similar to 
the no-overshoot models except that the convective core is larger 
(Figure 2). There is no mixing beyond the convection zone by de¬ 
sign. One difference is that the growth of the core can eradicate 
some of the remnants of the previous core-flash burning. Overall, 
the period spacing pattern is similar to the no-overshoot models 
(Figure 13). Importantly for the AIIi discrepancy, however, the 
larger convective core also increases the mean AP, while modes 
of the same radial order have a lower frequency. 


Figure 16. Pulsation properties of a 1 Mq classical semiconvection model 
with exactly V,.ad = Vad outside the convective core (black) and other¬ 
wise identical models where the composition discontinuity at the outer edge 
of the semiconvection region has been smoothed over 0.01 Mq (orange 
dashes) and 2 X 10“^ Mq (cyan). The region beyond the semiconvection 
zone has been homogenized in each model (see Section 2.2). The model 
without smoothing has AIIi = 273 s, the model with fine smoothing has 
AIIi = 271s, and the model with broad smoothing has AIIi = 264 s. 
These models have approximately R = 10.4 Rq, Teff = 4640 K, and 
t^max — 28 /iHz. 


3.5 Matching ensemble AIIi observations 

In this section we compare the inferred AIIi from the population of 
Kepler field stars with predictions from evolution models. We have 
chosen two representative masses: 1 Mq which experiences typical 
evolution for a red clump star (due to the uniformity of H-exhausted 
core mass at the flash), and 2.5 Mq which is massive enough to 
undergo core He-ignition in non-degenerate conditions (i.e. avoid 
the core flash) and then move to the so-called secondary clump in 
the HR diagram. We compare the models and observations with 


probability density functions P(Ani) constructed by the addition 
of Gaussian functions according to 


P(Ani) 


1 1 (AHi — Ani,i)^ 

i=l ^ 


( 20 ) 


where AHi,! represent each value from observations, or in the case 
of models, calculations at 1 Myr intervals, and N is the total of 
number of observations or calculated values. We use a standard de- 
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Figure 17. Comparison between the pulsation properties of a 1 Mq model 
with the classical semiconvection with exactly Vrad = Vad outside the 
convective core (black) and a model produced by the evolution code using 
our new mixing scheme (orange; described in Section 2.3.3). The region 
beyond the semiconvection zone has been homogenized in both models. 
The classical semiconvection model has AIIi = 273 s while the model 
from the evolution code has Afli = 268 s. The model from the evolution 
code has R = 11.1 R 0 , Tgff = 4600 K, and Umax = 31 //Hz (see Figure 16 
for the properties of the other model). 


viation of cr = 4 s and cr = 8 s for the 1 Mq and 2.5 Mq cases, 
respectively (Figures 18 and 19). 

In the 1 Mq case the maximal-overshoot models have the 
highest AIIi, followed by the standard-overshoot, semiconvec¬ 
tion, and no-overshoot models. The spreads of the ATEi probability 
density functions for the semiconvection and overshoot cases are 
smaller than is observed, and offset to lower values, as shown in 
Figure 18. In contrast, the spread for the maximal-overshoot 1 Mq 
model appears too broad, especially considering that we have com- 




Core helium burning age (Myr) 


Figure 18. Upper left panel: evolution of 1 Mq CHeB models with dif¬ 
ferent mixing schemes (no overshoot, standai'd overshoot, semiconvection, 
and maximal overshoot; in black, cyan, orange, and magenta respectively) 
in the Au — AIIi plane. Markers are at 10 Myr intervals. Determinations 
for Kepler field stars (grey dots) are from Mosser et al. (2014), and are lim¬ 
ited to those with reported mass 0.8 < M/Mq < 1.25. Upper right panel: 
probability density curves (Equation 20) for models in the upper left panel 
(same colours), standard overshoot with AIIi computed using only the re¬ 
gion outside the partially mixed zone (cyan dots), standard overshoot with 
increased Mhc (cyan dashes), and observations (grey dashes). Lower panel: 
surface luminosity evolution for the models in the upper left panel. 


puted single evolution sequences rather than a population which 
would widen the distribution. 

Every one of the four low-mass models appears to spend too 
much time with a low Alfi. Two possible resolutions are i) an in¬ 
creased H-exhausted core mass at the flash, which increases AIIi 
(dashed lines in Figures 8 and 18), or ii) that there is a difficulty in 
observationally determining AIIi for stars that have recently be¬ 
gun core helium burning (discussed in Section 3.3). Both of these 
affect the beginning of the CHeB, when AHi is lowest. The fact 
that this discrepancy exists even for the maximal-overshoot run, 
when the convective core is the largest possible, suggests that the 
treatment of convective boundaries cannot be the sole reason for it. 
In addition. Figure 19 shows that there is no evidence that this prob¬ 
lem exists for any of the higher-mass models. These more massive 
models do not experience the core flash, do not ascend the RGB 
to as high luminosity, and have more luminous hydrogen burning 
at the beginning of the CHeB phase, and thus would be unaffected 
by the proposed resolutions. In Figure 18 the appearance of the 
discrepancy at low AHi is worsened for the semiconvection and 
standard-overshoot runs by the slow decrease in AHi towards the 
end of CHeB. This alone cannot explain the discrepancy, however. 
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Figure 19. Upper left panel: evolution of 2.5 Mq CHeB models with dif¬ 
ferent mixing schemes (no overshoot, standai'd overshoot, semiconvection, 
and maximal overshoot; in black, cyan, orange, and magenta respectively) 
in the Av — AIIi plane. Markers are at lOMyr intervals. Determinations 
for Kepler field stars (grey dots) are from Mosser et al. (2014), and are lim¬ 
ited to those with reported mass 2 < AT/Mq < 3. Upper right panel: 
probability density curves (Equation 20) for models in the upper left panel 
(same colours), observations (grey dashes), and standard-overshoot model 
with Aril computed using only the region outside the partially mixed zone 
(cyan dots). Lower panel: surface luminosity evolution for the models in the 
upper left panel. 


because it is still present for sequences that do not undergo this 
slow drop in AIli late in CHeB (e.g. the dotted curve in Figure 18; 
discussed later in this section). 

The maximal-overshoot model is the only one of the four with 
different mixing prescriptions that can reach AHi values consistent 
with the hulk of the low-mass observations. Among the remaining 
cases, the standard-overshoot model is closest to the observations. 
The shape of its AHi probability density function also looks rea¬ 
sonable, except that it is offset by at least 25 s. Even a substantial 
increase in the H-exhausted core mass AMhc = 0.025 Mq (the 
most permitted by Catelan et al. 1996) at the start of core helium 
burning is not enough to match the entire observed AHi range. In 
that case it shifts the AHi probability density function higher by 
around 20 s. 

The 1 Mq semiconvection sequence has a lower AHi than our 
standard-overshoot case, by around 10 s. This is despite the similar 
evolution of i?cc and Mhc which strongly influence AHi (see Sec¬ 
tion 3.2). It is also evident from Figure 2 that the evolution of the 
size of the partially mixed region is similar for both sequences. We 
therefore attribute the difference in AHi to the way the composi¬ 


tion always varies smoothly in the semiconvection case, increasing 
N over a large interval in radius instead of over sharp spikes. 

We have performed an explicit test of the effect of the steep¬ 
ness of composition profiles on AHi. The three models in Fig¬ 
ure 16 are identical except for the composition near the edge of 
the semiconvection zone at r ~ 2.2 x 10® cm. The buoyancy fre¬ 
quency is nearly identical elsewhere in the structure (Figure 16b) 
so any difference in AHi must be due to the composition smooth¬ 
ing. In this case, smoothing the discontinuity over Am = 0.01 Mq 
decreases AHi by 9 s. 

This effect is also apparent in Figure 20, where smoothing the 
edge of the fully mixed core increases AHi. It can be seen in panel 
(a) that this smoothing increases the width of the peak in N (in 
the log scale) by more than it reduces its height. This increases the 
area under the curve, which reduces AHi. This is evident when we 
rewrite the integral in the asymptotic solution for AH^ in terms of 
log r to get 




-1 

r 1 

271^ 

/fdr 

27r2 

J A^d Inr 

y/t(l+l) 



( 21 ) 

The 1 Mq sequence without overshoot has the lowest AHi. 
AHi stays around 250 s, around 50 s below the bulk of the obser¬ 
vations, for the entire CHeB phase. This can be attributed to the 
lack of growth of the convective core (Figure 8c). The evolution 
of AHi in our 1 Mq and 2.5 Mq sequences without overshoot is 
almost identical to the corresponding models (also without over¬ 
shoot) from MESA (Stello et al. 2013). 

We have emulated the effect of mode trapping on the more eas¬ 
ily observable (non-trapped) modes in the standard-overshoot mod¬ 
els (discussed in Section 3.4.2) by excluding the partially mixed 
region from the calculation of AHi (dotted curve in Figure 18b). 
The impact of this is increasingly significant as core helium burn¬ 
ing progresses and the partially mixed region grows. This makes 
the AHi evolution very similar to that resulting from the maximal- 
overshoot scheme, except that it slightly exceeds the observed val¬ 
ues (by less than 10 s) near the end of core helium burning. How¬ 
ever, this is late in CHeB when this crude approximation of the 
effects of mode trapping is least valid, because the mode trap¬ 
ping cavity, and consequently the fraction of modes that become 
trapped, is large (making a neat fit in the period echelle diagram 
difficult; see Section 3.9). At the other extreme, AHi is still too 
low in the early stages of CHeB compared to the observations. 

The CHeB lifetime of the 1 Mq no-overshoot model is by far 
the shortest, followed by the maximal-overshoot case (Figure 18c). 
The semiconvection and overshoot sequences have nearly identi¬ 
cal lifetimes. The surface luminosity of the models is independent 
of the mixing scheme when they are still burning helium in the 
core (the variation in logL/LQ is less than 0.01). The relative 
energy generation rates from hydrogen and helium burning differ 
by a little more. The semiconvection model has more luminous H 
burning than the other sequences, while the no-overshoot model 
has the most luminous He burning, and the standard-overshoot and 
maximal-overshoot sequences are almost identical until the occur¬ 
rence of a core breathing pulse after 98 Myr. 

Star counts in globular clusters can be used as a constraint on 
the mixing scheme. This is because the CHeB lifetime is dependent 
on the amount of helium that is transported into the core. Models 
that consume less helium during CHeB have more helium that must 
be burnt during subsequent shell helium burning, and consequently 
have a longer early-AGB lifetime (compare the swift exit from the 
red clump of the no-overshoot sequence with its sluggish ascent of 
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the early-AGB in Figure 18c). The parameter R 2 = nAGn/ncHcB 
(the number ratio of observed AGB to CHeB stars) for globular 
clusters is thought to correspond to the ratio of the respective phase 
lifetimes. Caputo et al. (1989) argued that models with semicon¬ 
vection, but without breathing pulses, give the best fit to observa¬ 
tions of the globular cluster M5. In their models, the suppression 
of breathing pulses (by not allowing the growth in the convection 
zone if it would increase the central helium abundance) increased 
i ?2 from 0.10 to 0.14 or 0.15 (depending on the extent of core-flash 
burning), matching observations. 

We have computed R 2 for our models by considering lumi¬ 
nosity bins comparable to the observed range in metal-rich globular 
clusters. We have defined the CHeB lifetime to be when log L/L© 
is within 0.1 of its mean value before core helium depletion, and 
the AGB to be when logL/L 0 is no more than 1.0 higher than 
the CHeB range. We find values of R 2 = fAOB/fcHeB of 0.110, 
0.113, and 0.117 for the standard-overshoot, semiconvection, and 
maximal-overshoot schemes respectively. This would make them 
practically indistinguishable from one another by observations of 
star clusters. In contrast, the no-overshoot model has R 2 = 0.743, 
which is a difference that could easily be detected. We will ad¬ 
dress constraints from star counts in the next paper in this series 
(by computing less massive and more metal-poor models relevant 
to Galactic globular clusters; Harris 1996). 

There are a number of common trends between the 1 Mq (Fig¬ 
ure 18) and 2.5 Mq (Figure 19) models. In the 2.5 Mq runs the 
mixing scheme has a very similar effect on mean AHi, CHeB life¬ 
time, and the H- and He-buming luminosity. We also find a very 
similar effect from our emulation of mode trapping in the standard- 
overshoot model (dotted curve Figure 19). The probability density 
functions for the 2.5 Mq models are very similar in shape to the 
those for the 1 Mq models, except that they cover a larger range 
of AHi. The more substantial increase in AHi during their evo¬ 
lution can be explained by the greater extent of the growth of the 
H-exhausted core (roughly 0.2 Mq compared with 0.05 Mq for the 
1 Mq runs), the importance of which was shown in Section 3.2. 

Compared to the lower-mass case, the agreement with ob¬ 
servations is markedly better for the 2.5 Mq sequences, with the 
exception of the no-overshoot model. The semiconvection and 
standard-overshoot models, however, still do not reach the high¬ 
est AHi observations. In contrast, the AHi evolution for the mode 
trapping and maximal-overshoot sequences match each other even 
more closely, and both exceed the highest observed values by con¬ 
siderably more than does the 1 Mq mode trapping case. The com¬ 
parison between observations and models, however, is more com¬ 
plex than for the low-mass case. More of the increase in AHi is 
due to the growth of the H-exhausted core, and we are comparing 
the models to a population more diverse in mass and smaller in 
number. Therefore it would be imprudent to draw strong conclu¬ 
sions about the mixing from this sample. We note that our mod¬ 
els do not match the observed Au (but do match the shape of the 
population’s distribution in Av — AHi space). This is not prob¬ 
lematic because Av can easily be decreased by adjusting (in this 
case reducing) the MLT mixing length parameter, without affect¬ 
ing AHi. Finally, we note that by the end of core helium burn¬ 
ing, the 2.5 Mq models are considerably more luminous than at the 
beginning (by around a factor of 2 , apart from the shorter-lived no¬ 
overshoot model; Figure 19). This could introduce an observational 
bias for the secondary clump towards more luminous evolved stars, 
which have higher AHi (for all of the mixing schemes we have 
examined). Accounting for such a bias would help to resolve the 
excess of predicted low-AHi stars that is apparent in Figure 19b. 


In Section 3.4 we demonstrated how mode trapping may lead 
to an overestimation of AHi. This is made possible because only 
a subset of mixed modes can be detected. Here we briefly con¬ 
sider how mode trapping affects the period spacing between the 
pairs of modes that are most likely to be detected, i.e. those with 
low inertia. For this, we compare a standard-overshoot model with 
mode trapping (Figure 14) to a maximal-overshoot model without 
mode trapping (Figure 13). Although these two models have dif¬ 
ferent AHi (281 s and 314 s, respectively), the respective values 
determined from the period echelle diagram, APcch = 315 s (Fig¬ 
ure 15) and APcch = 316 s, are nearly identical. The average AP 
between all modes with 20 /rHz < 1 / < 40 /rHz for the standard- 
overshoot model is 270 s, which increases to 293 s if all of the (pre¬ 
sumably undetectable) trapped modes are excluded, compared with 
295 s for the maximal-overshoot model. When this calculation is re¬ 
stricted to the six pairs of modes closest to each low inertia trough 
(e.g., near v = 24 /rHz in Figure 14) we find AP = 275 s and 
AP = 277 s for the standard-overshoot and maximal-overshoot 
models, respectively. If we restrict the count to sets of four pairs of 
low-inertia modes instead of six we again find that the two mod¬ 
els have a similar average AP, except that it is reduced further, 
by 9 s in both cases. Moreover, the average frequency spacing be¬ 
tween these troughs is the same for both models. This indicates 
that knowing the typical observed AP would not assist with the 
detection of mode trapping. It also supports our suggestion in Sec¬ 
tion 3.4.2 that modes that are not trapped behave as though the 
buoyancy cavity is smaller than its true size, i.e. it excludes the 
semiconvection/partially mixed region with > 0 that is sur¬ 
rounded by a steep composition gradient that can trap modes. Be¬ 
cause the observationally determined AP depends on how many 
modes are detected it is difficult to compare these results to the 
average or median AP found in populations of CHeB stars (e.g. 
Mosser et al. 2011; Stello et al. 2013). 

3.6 The effect of the boundary of the convective core on 
pulsations 

In Figure 20 we examine the pulsations resulting from a structure 
with a fully mixed convective core but with a smooth composition 
profile at its boundary. This structure (which was produced by arti¬ 
ficially smoothing according to Equation 17) is interesting because 
of the physical implausibility of the core boundary in the maximal- 
and no-overshoot models (e.g. the magenta and black lines in Fig¬ 
ure 3b). In these models there is a true composition discontinuity 
where material that is strongly convectively unstable (Vrad Vad) 
does not partially mix with the material directly adjacent to it. 

All four models in Figure 20 with different composition pro¬ 
files have a comparable period spacing over much of the fre¬ 
quency range shown. The model in black with the largest partially 
mixed region (Am = 0.01 Mq), however, has several frequency 
ranges where the modes are more closely spaced (e.g., at around 
35 and 45 /rHz). This behaviour is also seen near zz = 38 ^Hz 
and V — 36 /jHz, for the models with Am = 0.002 Mq and 
Am = 5 X 10~^ Mq, in orange and cyan, respectively. 

These interruptions to the regular AP pattern shown in Fig¬ 
ure 20 c are more prevalent when the composition is smoother. 
This can be explained by the increasing buoyancy radius (see Sec¬ 
tion 3.2) of the mode trapping region enclosed by smoother com¬ 
position gradients (detailed in Section 3.5). This trapping region is 
adjacent to the convective core, however, so it always has a small 
buoyancy radius and therefore has little effect on mode inertia (Fig¬ 
ure 20d). This small buoyancy radius of the trapping region also ex- 
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Figure 20. Pulsation properties of 1 Mq models with a small partially 
mixed regions outside the convective core. In all four cases the composi¬ 
tion profile in the overshoot region has been set according to Equation 17 
with Am/M 0 values of 2 X 10“®, 5 X 10“^, 0.002, and 0.01 in ma¬ 
genta, cyan, orange, and black, respectively. In the same order AIIi is 
287 s, 294 s, 297 s, and 300 s. The composition between the convective core 
and H-buming shell has been homogenized. These models have approxi¬ 
mately R = 10.4 R 0 , Teff = 4770 K, and t/max = 29 /rHz. 


plains the long periodicity (in radial order) in its effect on AP, be¬ 
cause this gives a large An according to Equation 19. This contrasts 
with standard-overshoot and semiconvection models (e.g. Figures 
14 and 16), where a similar buoyancy feature is surrounded on both 
sides by a g-mode cavity which triggers mode trapping at regular 
intervals in radial order n with smaller An. Finally, we note that 
the disruption to the regular period spacing caused by a composi¬ 
tion gradient at the edge of the convective core appears most ob¬ 
vious for gravity-dominated modes, which are the most difficult to 
detect. 


3.7 Subdwarf B models 


We have also tested the effect of the core mixing scheme in subd¬ 
warf B models. At the beginning of core helium burning in these 
runs we homogenized the composition between the H-shell and 
the convective core to remove traces of core-flash burning. In each 
model we also set the helium mass fraction Y, at the H-exhausted 
core boundary according to 


AV 

Y (m) = Fsurf H —— 1 + cos 


/ m - Mhc 
V Am 



( 22 ) 


where AK is the difference between the surface and interior helium 
abundance, Mue is the mass of the H-exhausted core, m is the mass 
coordinate, and we have chosen Am = 0.002 M©. We chose this 
smooth profile because we are only interested in the effect of the 
composition profile at the boundary of the convective core. We set 
the total mass and mass of the shell to match the mass, gravity, and 
effective temperature typical of the stars in the Reed et al. (2011) 
sample. This ad hoc approach is obviously inadequate for precision 
studies of particular stars (e.g., those found in Charpinet et al. 2011; 
Van Grootel et al. 2013b), but suits our purpose here. 

The results of the pulsation calculations for the four models 
with different mixing schemes are presented in Figure 21. The ap¬ 
pearance of the AP pattern for each model is broadly similar, ex¬ 
cept for a few trapped modes in the standard-overshoot model (Fig¬ 
ure 21). These modes have much higher inertia than their neigh¬ 
bours and are more closely spaced in period. This behaviour is 
similar to our red clump standard-overshoot model (Figure 14) and 
the semiconvection model that includes a region with a stabilizing 
molecular weight gradient at the edge of the semiconvection zone 
(Figure 16). 

There is a substantial difference in the mean AP between the 
four different mixing cases. AP spans a range of around 60 s, with 
the no-overshoot sequence having the lowest average value, fol¬ 
lowed by the semiconvection case, then the maximal-overshoot and 
standard-overshoot cases (if the trapped-modes are excluded). Like 
its more massive counterparts, the sequence without overshooting 
has a lower AHi than is observed. In this case it is more than 10 s 
too low to match any of the observations reported by Reed et al. 
(2011), which is an especially strong constraint because the model 
in Figure 21 is from the stage of CHeB when AHi is near its max¬ 
imum. The range of AP between pairs of modes found for KIC 
5807616 spans less than 30 s (Reed et al. 2011) which is consistent 
with the range for our models (except near the high-inertia modes 
in the standard-overshoot case). 

Recently, 0stensen et al. (2014) found evidence for mode 
trapping in KIC 10553698A, an sdB star in the Kepler field, by clas¬ 
sifying 1 = 1 and £ = 2 modes. They identified both the C-O/He 
and He/H transition zones as possible origins of the mode trapping 
and highlighted the resemblance of the period spacing pattern to ex¬ 
isting theoretical calculations, e.g. Fig. 3 in Charpinet et al. (2002). 
Due to the use ofq = log [1 — m/M] for the horizontal axis of that 
figure, the structure near the core is difficult to discern, but it ap¬ 
pears that in their “evolutionary model” there is a relatively smooth 
buoyancy peak near where partial mixing can occur in our models. 
In our standard-overshoot model the mode trapping is certainly a 
result of the sharp composition gradient at the edge of the partially 
mixed zone outside the convective core. Moreover, that model’s 
pulsations bear perhaps an even more remarkable similarity to the 
observations shown in grey in Figure 21 (keeping in mind we made 
no attempt to match the frequencies). The trapped modes in our 
model, however, reside deep within the core, so their observability 


MNRAS 000, 000-000 (0000) 





















Astewseisrnology of core helium burning stars 19 


Table 1. Properties of the evolution sequences. The models have standard 
overshoot (SO), no overshoot (NO), semiconvection (SC), and maximal 
overshoot (MO). The triple-a and ^^C(o, 7)^®0 reaction rates are denoted 
by fca ™d rsQ,, respectively. AIIi mean is the average AIIi value in the 
CHeB phase which has duration rcHeB- The initial and final H-exhausted 
core masses for the CHeB phase are denoted by Mhcj and Mne.f respec¬ 
tively. 


Model notes 

M 

(Mq) 

Mixing AHl^mean rcHeB ^He.f 

(s) (Myr) (Mq) (Mq) 


1 

SO 

267 

135.0 

0.466 

0.503 


1 

SC 

258 

129.6 

0.467 

0.503 


1 

NO 

204 

79.5 

0.467 

0.499 


1 

MO 

293 

119.7 

0.467 

0.499 


2.5 

SO 

232 

221.1 

0.331 

0.509 


2.5 

SC 

225 

227.5 

0.331 

0.514 


2.5 

NO 

182 

170.0 

0.331 

0.475 


2.5 

MO 

251 

216.8 

0.331 

0.502 

AMne = 0.025 Mq 

1 

SO 

285 

99.8 

0.493 

0.520 

AMne = 0.025 Mq 

1 

MO 

304 

100.4 

0.493 

0.519 

rca, r3a X 2 

1 

SO 

276 

137.5 

0.462 

0.494 

rca X 2 

1 

SO 

275 

132.8 

0.466 

0.503 

rsa X 2 

1 

SO 

270 

119.9 

0.462 

0.491 

y-1-0.1 

1 

SO 

272 

112.0 

0.448 

0.548 

[Fe/H] - 1.0 

1 

SO 

285 

112.4 

0.467 

0.523 


is uncertain. Finally, the theoretical AP for our non-trapped modes 
is nearly an exact match for the AP between most observed modes 
in KIC 10553698A, which suggests that the size of the convective 
core in the standard-overshoot model is reasonable. 

We have also examined the effect of core flash phase burn¬ 
ing on the pulsations in our sdB models. In Figure 22 we demon¬ 
strate that this effect is strongly dependent on the smoothness of 
the remaining composition profile. The four models that we use to 
test this include a model without the discontinuity from the core 
flash (constant composition), and others with sine wave composi¬ 
tion profiles (Equation 17), with Am set as 0.01 Mq, 0.001 M 0 , 
and 10“® Mq. The consequences for the computed frequencies are 
increasingly apparent for models with sharper composition pro¬ 
files. The mode period spacing and inertia for the model with 
the smoothest composition profile (spread over 0.01 Mq) is nearly 
identical to the model with a constant composition. By comparison, 
the model with the chemical profile spread over 0.001 Mq shows 
up to four times the period spacing variation for high frequency 
modes {P < 1.3 x 10“^ s) and mode to mode variation of almost 
100 s at lower frequency (where the two smoother models show al¬ 
most constant AP). The model with the sharpest composition pro¬ 
file shows a pattern similar to the Am = 0.001 Mq case, except 
with more extreme variation in AP. In both cases the amplitude of 
this variation oscillates, with a period of around 1.8 x 10^ s. We 
also note the similarity between the period spacing pattern of these 
two models and the model by Charpinet et al. (2014). 

Overall, these results suggest that the possibility of using pul¬ 
sations to determine whether a low-mass CHeB star has experi¬ 
enced the core flash depends principally on how discontinuous is 
the composition profile it has left behind. 

3.8 Dependence on input physics 

In Section 3.2 we showed that 1 Mq standard-overshoot models 
need an increase in the H-exhausted core mass at the flash of more 
than AMne > 0.025 Mq to match the range of AHi reported for 




Figure 21. Seismic properties of synthetic subdwarf B models and an ob¬ 
served pulsator in the Kepler field. Upper panels: Brunt-Vaisala frequency 
N and helium mass fraction Y for models with different mixing prescrip¬ 
tions. Lower panels: £ = 1 mode spacing AP and inertia. The observations 
are of £ = 1 modes classified by 0stensen et al. (2014) for KIC 10553698A 
(thick grey lines and squares). The models were generated according to the 
method outlined in Section 3.7. They have M = 0.475 Mq, solar metallic- 
ity, and yi;ent = 0.4. The models have no overshoot (black), standard over¬ 
shoot (cyan), semiconvection (orange), and maximal overshoot (magenta). 
These models have Afli of 222 s, 245 s, 238 s, and 269 s, respectively, and 
approximately R = 0.20 Rq and T^ff = 27000 K. 


low-mass CHeB stars by Mosser et al. (2014). The effect of un¬ 
certainties in the input physics on Mhc at the core flash has been 
examined in detail previously (e.g. Catelan et al. 1996). Some of 
these uncertainties are not important to subsequent CHeB evolu¬ 
tion. For instance, the expansion of the core during the flash phase 
decreases both its rotation rate and neutrino emission so these ef¬ 
fects need only be considered in light of how they affect the core 
mass at the flash. In contrast, helium burning reaction rates and 
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Figure 22. Comparison of the pulsation properties of sdB models with dif¬ 
ferent composition discontinuities left behind by the core flash. One model 
(orange) has a constant composition between the convective core and the 
H-buming shell. The models in black, magenta, and cyan have chemi¬ 
cal gradients over 0.01 Mq, 0.001 Mq, and 10“® Mq, respectively. Each 
composition profile is set according to Equation 17. These models have 
R = 0.20 Rq and Teff = 27400 K. 


initial composition also affect the later evolution, including AIIi 
(Table 1). 

Doubling the ^^C(a,7)^®0 reaction rate increases the aver¬ 
age Afli during CHeB by 8 s. Once there is enough carbon in 
the core (Ac ^ 0.1) the ^^C(a, 7)^®0 reaction proceeds more effi¬ 
ciently, slowing the rate of increase of the central temperature and 
density (and therefore also the triple-Q rate). This increases the 
convective core radius and consequently Afli, but does not signifi¬ 
cantly affect the CHeB lifetime. In contrast, increasing the triple-a 
rate reduces the core temperature and density from the beginning 
of CHeB. Although this tends to increase AHi, it is offset by the 
lower Mhc, which starts smaller and grows more slowly due to the 


consequently reduced hydrogen-burning luminosity. This results in 
only a 3 s increase in the average AHi during the CHeB phase. In 
the relevant conditions the uncertainty in the triple-Q rate is less 
than 15 per cent while for ^^C(q, 7)^®0 rate it is around 40 per 
cent (Angulo et al. 1999), and more recent data favours the lower 
limit (Xu et al. 2013). Taking both of these uncertainties into ac¬ 
count, they could together only account for around a 5 s change in 
AHi, considerably less than the size of the disparity between stan¬ 
dard models and observations (of around 30 s). 

We have also tested the consequences of varying the initial 
composition. Increasing helium raises the average AHi, but the 
dependence is weak: a large increase of AY = 0.1 only increases 
the average AHi by 5 s. This may be attributed to the more rapid 
growth of the H-exhausted core during CHeB compared to the stan¬ 
dard case, making it 0.04 Mq larger at core helium exhaustion. This 
is partly offset, however, by the lower H-exhausted core mass at he¬ 
lium ignition, limiting the increase in average AHi. Reducing the 
metallicity by a factor of 10 increases the average AHi during the 
CHeB phase by 18 s. This is due to a reduction in the heavy ele¬ 
ment opacity (which we confirmed by evolving an [Fe/H] = — 1 
model but with solar heavy element opacity; which had a negli¬ 
gible effect on AHi). This initially increases the helium burning 
rate and consequently AHi. The hydrogen burning rate increases 
even more substantially, which further increases AHi by acceler¬ 
ating the growth of the H-exhausted core. Composition, however, 
is not likely to be the cause of the AHi discrepancy because the 
stars in the Mosser et al. (2014) sample are typically around solar 
metallicity (Pinsonneault et al. 2014), consistent with the models in 
Section 3.5. Indeed, none of these factors, nor any reasonable com¬ 
bination of them, can explain why the 1 Mq standard-overshoot run 
fails to match the observations. 

3.9 Late-CHeB and early-AGB models 

A number of authors have identified possible late-CHeB and AGB 
stars in the Kepler field through seismology. Mosser et al. (2012) 
found five stars with the same Av as the low-mass CHeB group 
but with lower AHi (around 250 s) and posited that these stars 
have exhausted helium in their cores. Corsaro et al. (2012) identi¬ 
fied several members of the open clusters in the Kepler field (NGC 
6811, NGC 6819 and NGC 6791) that are likely to be evolved red 
clump stars because they have similar 1=1 AP to the majority of 
clump stars, but have lower Av. In their examination of field stars, 
Kallinger et al. (2012) suggested that lower-Ai/ stars belong to the 
early-AGB. These stars also have a distinct central radial {I = 0) 
mode phase shift, which can be attributed to a difference in the 
structure of the convective envelope (Christensen-Dalsgaard et al. 
2014). Although it is not examined in this paper, it would be in¬ 
teresting to determine if and how this phase shift depends on the 
CHeB mixing scheme. 

Our models disagree with the earlier suggestion by 
Mosser et al. (2012) that red clump stars with a low AHi but typ¬ 
ical Av can be explained as being post-CHeB (they are now clas¬ 
sified with the other red clump stars in Mosser et al. 2014). These 
are unlikely to be post-CHeB because every one of our low-mass 
models - irrespective of mixing scheme - shows a decrease in Av 
when AHi begins to decrease, which occurs prior to central he¬ 
lium exhaustion (Figure 8). This causes them to move away from 
the location of the suspected post-CHeB stars in Av — AHi space 
(Figure 18a), which is in agreement with the MESA models without 
overshoot shown in Figure 4b in Stello et al. (2013). 

In Figure 23 we show a standard-overshoot model before 
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and after core helium exhaustion, separated by 160 kyr. During 
this period there is a rapid increase in luminosity (logL/L© in¬ 
creases from 2.029 to 2.117), a decrease in l\v (from 1.80/rHz 
to 1.47 fjHz), and a decrease in AIIi (from 153 s to 99 s; it then 
drops to 65 s after a further 1 Myr). This sudden decrease in pe¬ 
riod spacing has also been shown for models computed with MESA 
(Stello et al. 2013). If the region enclosed by the outer edge of the 
partially mixed zone (dashed lines in Figure 23) is excluded from 
the calculation of ADi (to emulate the effect of mode trapping) the 
drop in period spacing is less severe (from 251 s to 234 s). This, 
however, still suggests that the high AP (~ 250 s) stars, identified 
as possible members of the AGB by Kallinger et al. (2012), are in 
fact still CFIeB stars. Similarly, the relatively low Au open cluster 
stars identified as evolved red clump stars in Corsaro et al. (2012) 
appear to be correctly classified, while the one suggested early- 
AGB star in NGC 6811 is probably also in the late-CHeB phase. 
Lastly, the position in Au — AIIi space of our late-CHeB models, 
and those from Stello et al. (2013), is generally consistent with the 
observed group marked by ‘A’ in Fig. 1 in Mosser et al. (2014). We 
suggest that care should be taken when describing these stars, be¬ 
cause “red clump” and “core helium burning” are not interchange¬ 
able terms. According to models, stars leave the red clump when 
they are still burning helium in the core. 

Determinations of AIIi from observations of stars near core 
helium exhaustion could be very uncertain if there is mode trap¬ 
ping in the partially mixed region. This is because late in CHeB the 
buoyancy radius of the partially mixed region, where the modes are 
trapped, becomes large compared to the total buoyancy radius. In 
the pre- and post-core helium exhaustion models in Figure 23 the 
partially mixed regions account for 41 per cent and 68 per cent of 
the total buoyancy radius respectively. There are thus few modes of 
consecutive radial order that are both not trapped, unlike the model 
in Figure 15 for instance. The extensive mode trapping in these 
models would make it difficult to accurately determine AIIi from 
observations, but also make it unlikely they could be interpreted as 
having an erroneously high AIIi from the period echelle diagram. 

Core breathing pulses (GBP) only occur in the standard- 
overshoot model (note the rapid increases in central helium abun¬ 
dance that begin after 98 Myr in Figure 8). GBP do, however, occur 
in each of the remaining models if the mixing scheme is changed 
to standard overshoot late in CHeB (when the central helium abun¬ 
dance is y = 0.1). This demonstrates that CBP are prevented by 
the mechanics of each mixing prescription rather than by the very 
different late-CHeB structures they eventually produce. An exam¬ 
ple of the divergence of the internal composition is shown in Fig¬ 
ure 24c and is discussed below. Although these structural differ¬ 
ences do not prevent CBP, they do affect the magnitude of them; a 
larger convective core, or the existence of a partially mixed region 
outside it, reduces the amount of helium transported into the core 
by the breathing pulses. 

Finally, in Figure 24 we show four examples of late-CHeB 
models that were evolved with different mixing schemes until they 
have central helium abundance of around Y = 0.01. The con¬ 
traction of the fully convective core is evident in ah but the no¬ 
overshoot model (note the convectively stable region near r = 
1.6 X 10® cm in the maximal-overshoot model in magenta). By 
this stage the partially mixed regions in the semiconvection and 
standard-overshoot models extend well beyond the edge of the 
maximal-overshoot core, whereas earlier in the evolution their sizes 
are comparable (e.g. Figure 3a). It is also clear that by the end of 
core helium burning those models have burned more helium than 
the maximal-overshoot case. By the end of CHeB the internal struc- 



r (lO^cm) 


Figure 23. The Brunt-Vaisala frequency N (upper panel) and helium abun¬ 
dance Y (lower panel) for a 1 Mq solai'-metallicity run with standard over¬ 
shoot before (black) and after (cyan) core helium exhaustion. The first 
model has central helium abundance ILcnt = 4 X 10““*. Between the two 
models 160 kyr elapses. In this time AIIi decreases from 153 s to 99 s, Au 
decreases from 1.80/rHz to 1.47/rHz, radius increases from IS.ORq to 
21.0 R 0 , Teff decreases from 4300 K to 4250 K, and i/max decreases from 
9 ^Hz, to 7 ^Hz. 

tures have diverged significantly enough to suggest that i) the mix¬ 
ing scheme could affect the early-AGB evolution, perhaps to an 
extent that is detectable in a large enough homogeneous population 
(e.g. globular clusters; which will be explored in a forthcoming pa¬ 
per) and ii) asteroseismic studies of the population of late-CHeB 
(such as that found in Mosser et al. 2014) and early-AGB stars may 
provide vital clues about CHeB evolution. 


4 SUMMARY AND CONCLUSIONS 

The asteroseismic detection of mixed modes in core helium burn¬ 
ing stars in the Kepler field offers an unprecedented insight into the 
internal structure of these stars. With the aim of better constrain¬ 
ing the models, we have investigated two discrepancies between 
the predicted asymptotic g-mode £ = 1 period spacing AHi from 
standard low-mass (1 Mq) stellar models and those reported for the 
Kepler field stars (Mosser et al. 2012, 2014): 

• The average value of AHi predicted by the models is signifi¬ 
cantly below the average inferred from observations (by more than 
25 s) and the models never reach the highest observed values of 
AHi. 

• The models spend more time with low values of AHi during 
core helium burning than is implied by the observed population 
(Figure 1). 

One possible source of these discrepancies could be that there 
are systematic problems with the internal stellar structure of stan¬ 
dard models. Indeed, it is well known, yet often ignored, that the 
models of this phase are uncertain (e.g. Figure 2). To explore these 
uncertainties in the light of the new asteroseismic observations we 
computed non-radial adiabatic pulsations and AHi for a suite of 
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Figure 24. Compaiison of the structure of four late-CHeB models with dif¬ 
ferent treatments of convective boundaries, each with ilent ~ 0.01. The 
mass range in the lower panel approximately corresponds to the radius 
range in the upper two panels. The colours are the same as in Figure 3. The 
standard-overshoot, semiconvection, no-overshoot, and maximal-overshoot 
models have AIIi of 208 s, 298 s, 195 s, and 247 s, respectively. These 
models have approximately R = 15 R©, Teff = 4440 K, and i/max = 
14 fiUz. 


core helium burning models with varying physical inputs and mix¬ 
ing algorithms. The stellar models were calculated with: 

(i) Four different mixing schemes (Section 2.3). 

(ii) Three different initial chemical compositions (Section 3.8). 

(iii) Altered He-buming reaction rates (Table 1). 

We found that varying the stellar composition or altering the 
He-buming reaction rates cannot reconcile the models and obser¬ 
vations (Section 3.8; Table 1). Three of the four mixing schemes 
also failed to increase AHi by the magnitude required. The only 
models that can match the average observed AHi values reported 
are those with large convective cores, such as those calculated with 
our newly proposed “maximal-overshoot” scheme (Section 2.3.4). 
In this scheme the extent of convective overshoot is adjusted so that 
it produces the most massive convective core possible. This treat¬ 
ment was implemented, however, only as a demonstration of the 
effect of a large convective core: we have not proposed any physi¬ 
cal basis for it. In the case of more massive stars (M > 2 M©), the 
smaller number of observations, and the fact that their H-exhausted 
cores grow substantially during the core helium burning phase, al¬ 
lowed us to only rule out the no-overshoot model (Figure 19). 

Another possible source of the AHi discrepancies is that the 
observations may be biased in some way. By comparing our non- 
radial adiabatic pulsation calculations against AHi across the suite 


of models, we identified a potential difficulty in inferring AHi 
from observations: any mode trapping that results from a con¬ 
vective region between two radiative zones (e.g. Section 3.2; Fig¬ 
ure 7), or a steep composition gradient at the outer boundary of a 
semiconvection or partially mixed region (e.g. Section 3.4.2; Fig¬ 
ure 15), increases the period spacing between most pairs of modes 
of consecutive radial order, and therefore the observationally in¬ 
ferred value of AHi. The difference between these values could 
explain much of the disagreement between standard models and 
observations (dotted curve in Figure 18). 

However, even after accounting for these two proposed reso¬ 
lutions to the discrepancy in average AHi values, the models still 
predict more core helium burning stars with low AHi (< 270 s) 
than observed. We suggested two possible remedies for this prob¬ 
lem: i) there may be a difficulty in observationally determining 
AHi for early core helium burning stars (when AHi is lowest) be¬ 
cause the sharp composition profile at the hydrogen burning shell 
causes the period spacing pattern to be highly irregular compared 
to more evolved models (Section 3.3), or ii) the mass of the he¬ 
lium core at the flash may be higher than predicted by standard 
models, thereby raising the initial AHi (Figure 18). Further in¬ 
formation about the selection effects in asteroseismic population 
studies, which are alluded to by Mosser et al. (2014), would help to 
establish the validity of the first point. At present, the possibility of 
unknown systematic biases in the observations limits our ability to 
use them to make firm conclusions about stellar evolution theory. 

We also investigated the dependence of AHi on the radius of 
the convective core, as shown by Montalban et al. (2013). We found 
that the relationship to the observed period spacing is more com¬ 
plicated than a simple linear relationship in a number of respects 
(Section 3.2). Furthermore, AHi is also dependent on the steep¬ 
ness of any chemical profiles outside the convective core, such as 
those found in the semiconvection zone (Section 3.5). 

The structure of low-mass CHeB stars is further complicated 
by the stabilizing chemical gradient left behind by helium burning 
during the core-flash phase. This can have a significant effect on 
the period spacing pattern, depending on the steepness of the gra¬ 
dient (Figure 12). In fact, in models taken directly from our evolu¬ 
tion code the mode trapping from the main discontinuity produced 
in the core-flash phase is the most important feature in the period 
spacing (e.g. Figure 13). This contrasts to composition gradients 
that may be created by overshooting (a small distance) from the 
convective core which have a more subtle effect on the period spac¬ 
ing (Figure 20). 

We also tested low-mass models that imitate sdB stars (Sec¬ 
tion 3.7). In these models we find the same dependence of AHi and 
AP from pulsation calculations on mixing scheme as for our solar- 
mass models. We also found that it may be difficult to use aster- 
oseismology to constrain sdB formation channels. This is because 
the effect of the composition discontinuity resulting from core-flash 
burning is smallest for the low radial-order modes (Figure 22) that 
are typically detected. We noted, however, that there may be other 
evidence from asteroseismology such as differences in the H-He 
transition region (Hu et al. 2008). 

In Section 3.9 we showed that our models of core helium ex¬ 
haustion suggest that early-AGB stars will not be found near the 
bulk of core helium burning stars in the Au - AHi diagram, inde¬ 
pendent of mixing scheme. This is because both AHi and Au have 
decreased by the time core helium burning ceases. This expands on 
the earlier finding by Stello et al. (2013) for models without con¬ 
vective overshoot. 
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Finally, although we have highlighted some possible expla¬ 
nations for the discrepancies in AIli, further work is needed to 
pinpoint the cause(s). In order to better gauge the extent of the 
problem, and therefore the merit of our proposed solutions, it is 
necessary to account for any selection bias in the observations. Our 
possible solution involving the mode trapping phenomenon affect¬ 
ing the observationally inferred AIIi values (e.g. Section 3.4.2) 
could be investigated by comparing models to specific frequency 
patterns observed. Constraints on the core mass at the helium flash 
and the mixing during the CHeB phase could be investigated by 
using the latest photometry of globular clusters - this is the subject 
of the next paper in this series. 
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